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Abstract Bats have been shown to serve as natural res- 
ervoirs for numerous emerging viruses including severe 
acute respiratory syndrome coronavirus (SARS-CoV). In 
the present study, we report the discovery of bat CoV genes 
in Indonesian Moluccan naked-backed fruit bats (Dobsonia 
moluccensis). A partial RNA-dependent RNA polymerase 
gene sequence was detected in feces and tissues samples 
from the fruit bats, and the region between the RdRp and 
helicase genes could also be amplified from fecal samples. 
Phylogenetic analysis suggested that these bat CoVs are 
related to members of the genus Betacoronavirus. 
Keywords Bat coronavirus - Fruit bats - Indonesia 
Coronaviruses (CoVs) are known to cause infections in 
both humans and animals, with varying symptoms ranging 
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from mild respiratory illness to severe infections resulting 
in death [1]. In addition to the respiratory tract, coronavi- 
ruses can also affect other organs in the body, including the 
gastrointestinal tract, liver, kidney, and brain of both 
humans and animals [2]. The pandemic of severe acute 
respiratory syndrome (SARS) in 2002 -2003 [3, 4] and the 
emergence of Middle East respiratory syndrome (MERS) 
[5] in 2012 are examples of human infections resulting in 
significant mortality. In the case of SARS, the mortality 
rate reached 10 % as a result of respiratory failure [6—8]. 

Coronaviruses are the largest known RNA viruses, with 
genomes up to 32 kb. They belong to order Nidovirales, 
family Coronaviridae. These viruses possess an envelope 
and a single-stranded, positive-sense RNA genome [9]. 
Currently, four genera have been described: Alphacoro- 


navirus, Betacoronavirus, Gammacoronavirus, and 
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Deltacoronavirus [10]. Members of the genera Alphacor- 
onavirus and Betacoronavirus are known to cause human 
disease, whereas those of the genera Gammacoronavirus 
and Deltacoronavirus are causative agents of animal dis- 
ease [11]. Bats are known to be reservoirs for many 
emerging infectious diseases that have zoonotic potential 
[12] and are believed to be the original source of SARS- 
CoV [13] and MERS-CoV [14]. In Southeast Asia, bats 
harboring coronaviruses have been discovered in a number 
of countries, including the Philippines [15, 16] and Thai- 
land [17]. However, there have been no studies focused on 
CoVs in Indonesia. Here, we report the discovery of bat 
CoVs in Moluccan naked-backed fruit bats, Dobsonia 
moluccensis, in Indonesia. 

We collected fruit bat samples from three regions in 
Indonesia: Paguyaman District, Gorontalo Province (in 
2012 [n = 37] and 2013 [n = 15]); Surabaya District, East 
Java Province (in 2012 [n = 3]); and Yogyakarta District, 
Yogyakarta Special Province (in 2012 [n = 19]). All ani- 
mal experiments were conducted in accordance with the 
ethical guidelines of the Animal Care and Use Committee 
of the Veterinary Teaching Hospital, Bogor Agricultural 
University. Fecal, tracheal swabs, and tissues samples 
including brain, lung, liver, spleen, and kidney from each 
bat were collected and placed into RNAlater (Life Tech- 
nologies, Carlsbad, CA). All samples were exported to 
Japan with the permission of the Directorate General of 
Livestock and Animal Health Services, the Ministry of 
Agriculture, Republic of Indonesia. 

Morphological features and nucleotide sequence analysis 
of both mitochondrial 16S rRNA and cytochrome b were 
conducted as described previously [18] to identify the 
species of each fruit bat. Feces and tissue samples were 
subjected to RNA extraction using TRIzol Reagent (Life 
Technologies) according to the manufacturer’s instructions. 
The RNAs were then screened for genomes of CoVs using a 
nested RT-PCR with primer sets that specifically amplify 


ORF1b 





Fig. 1 Schematic diagram of genes that were amplified for detection 
of CoVs from fruit bats based on coronavirus replicase gene. ORFla 
and ORF 1b are constituents of the replicase gene. ORFla, open 
reading frame la; RdRp, RNA-dependent RNA polymerase; Hel, 
helicase; ORF1b, open reading frame 1b; A, partial RdRp gene 
(455 bp); B, partial helicase gene (355 bp); C, region between 
amplified partial RdRp and partial helicase gene (1,335 bp); D, partial 
RdRp-helicase gene (2,145 bp). Arrows indicate the orientation of the 
primers used for RT-PCR 
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the conserved region of RdRp, as described previously [19] 
(region A in Fig. 1). The helicase region was amplified 
using degenerate primers specific for members of the genus 
Betacoronavirus [5] (region B in Fig. 1). PCR products of 
both the RdRp and helicase regions were gel-purified, and 
nucleotide sequencing was performed using a BigDye 
Terminator v.1.1/3.1 Cycle Sequencing Kit (Life Technol- 
ogies) and analyzed using an ABI PRISM 3130 Genetic 
Analyzer (Life Technologies). One-step RT-PCR was per- 
formed using a new primer set and a PrimeScript II High 
Fidelity OneStep RT-PCR Kit (TaKaRa, Tokyo, Japan) to 
amplify the region not covered in the first and second sets of 
PCR experiments (region C in Fig. 1). Based on alignment 
of the sequences obtained in this study, the forward primer 
5'-CACGCAACTTGTTGTAATGCGTCAGAGA-3’ and 
the reverse primer 5'-CACGTGCTTTTGCAGGCACTAT 
ACGAC-3’, corresponding to the RdRp gene at nucleotide 
position 14850-14877 and to the helicase gene at nucleotide 
position 16714-16740, respectively, of bat coronavirus 
(BatCoV) HKU9 (NC009021) were used. The fragments 
obtained, covering part of the RdRp and helicase genes 
(region D in Fig. 1), were sequenced. 

PCR-amplified nucleotide sequences were aligned and 
compared to the sequences of known CoVs using the 
MEGA 6 program. Nucleotide sequence identity values 
were calculated using GENETYX software ver. 10 
(GENETYX, Tokyo, Japan). Phylogenetic trees were 
constructed based on nucleotide sequences, using the 
maximum-likelihood method with the Tamura-Nei model 
and 1000 bootstrap replicates [20, 21]. 

Virus isolation was attempted by inoculation of African 
green monkey kidney (Vero E6), Yaeyama flying fox 
kidney (FBKT), and Leschenault’s rousette kidney (Dem- 
KT1) cells. Vero E6 cells were obtained from Prof. Ayato 
Takada (Research Center for Zoonosis Control, Hokkaido 
University, Hokkaido, Japan) [22], and FBKT and Dem- 
KT1 cells were obtained from Prof. Ken Maeda (Yamag- 
uchi University, Yamaguchi, Japan) [23]. Fecal and tissue 
homogenates that were positive for CoV RNA were used 
for inoculation of cells after clearance by centrifugation at 
1000xg for 5 min. Cells were incubated at 37 °C for 
30 min with fecal samples and for 1h with tissue 
homogenates. Inocula were replaced by fresh culture 
medium supplemented with 2 % FBS, 2 % antibiotic- 
antimycotic solution containing penicillin (10,000 U/ml), 
streptomycin (10 ug/ml) and amphotericin B (25 ug/ml) 
(Sigma, St. Louis, MO), and gentamycin (25 ug/ml). The 
supernatants of cultured cells were transferred to fresh cell 
monolayers each week. RNAs were extracted from culture 
supernatants and subjected to nested RT-PCR targeting the 
RdRp gene up until the fifth passage. The presence of 
cytopathic effect (CPE) was also checked daily until the 
end of the experiment. 


Coronavirus genomes in fruit bats in Indonesia 


The species of the captured bats were shown to be Dob- 
sonia moluccensis, Acerodon celebensis, Pteropus sp. 
(genetically close to Pteropus hypomelanus), and Pteropus 
vampyrus (Table 1). Nested RT-PCR with three out of a 
total of 74 fecal samples from Dobsonia species produced a 
fragment of 455 bp in length after removal of the primer 
region, which was expected size of the RdRp fragment 
(region A in Fig. 1 and Table 1). The samples were 
IFB2012-8F, IFB2012-13F, and IFB2012-17F. The RdRp 
gene is the one most commonly used for PCR amplification 
in CoV surveillance, since it contains the most regions that 
are conserved in all CoVs [24]. Sequence analysis demon- 
strated that all samples contained the sequence of the RdRp- 
gene-conserved motifs A (DYPKCD) and C (XSDD), which 
form the polymerase catalytic active site [24]. The partial 
RdRp sequences of IFB2012-13F and IFB2012-17F were 
identical. There were two nucleotide substitutions when the 
RdRp sequences of IFB2012-13F and IFB2012-17F were 
compared with the RdRp sequence of IFB2012-8F. 

In addition to the partial RdRp gene, a portion of the 
helicase gene was also amplified, producing a fragment of 
355 bp (region B in Fig. 1), from the three samples that 
were positive for PCR targeting the RdRp region. In the 
helicase gene, no nucleotide differences were found 
between IFB2012-13F and IFB2012-17F. However, two 
nucleotide substitutions were observed between IFB2012- 
8F and IFB2012-13F or IFB2012-17F. 

After sequencing the partial RdRp and partial helicase 
regions, we designed new primers to amplify the sequence 
between these regions (region C in Fig. 1). A 1335-bp 
fragment was detected in two samples (IFB2012-8F and 
IFB2012-17F). The sequence of a fragment (region D in 
Fig. 1), referred to as “partial RdRp-helicase”, was then 
constructed by incorporating overlapping sequences of par- 
tial RdRp, partial helicase, and the region in between. The 
nucleotide sequence comparison of partial RdRp-helicase 
from IFB2012-8F and IFB2012-17F showed 13 nucleotide 
substitutions. BLAST searches showed that IFB2012-8F and 
IFB2012-17F had the highest nucleotide sequence identity to 
BatCoV HKU9-10-2 and BatCoV KY06, at 86 % and 85 %, 


Table 1 Detection of 


Year Bat species 
coronaviruses in fruit bats by 
nested RT-PCR targeting RdRp 
gene in fecal samples 2012 Dobsonia mollucensis 


Acerodon celebensis 


Pteropus sp. 


Pteropus vampyrus 
Pteropus vampyrus 
2013 Acerodon celebensis 


Pteropus vampyrus 


Total 


respectively. Following removal of primer sequences from 
the partial sequence of the RdRp-helicase gene at the 5’ and 
3’ end, the length of the sequence was 2,116 bp, and a phy- 
logenetic tree comparing the sequence of the samples with 
known CoVs, including SARS-CoV and MERS-CoV, was 
constructed (Fig. 2). The phylogenetic tree showed that 
IFB2012-8F and IFB2012-17F formed a distinct branch that 
was closely related to BatCoV HKU9 (NC009021), HKU9-2 
(EF605514), HKU9-5-2 (HM211099), and HKU9-10-2 
(HM211101) from China [25] as well as BatCoV K Y06 from 
Kenya (HQ728483) [26]. All of these CoVs belong to the 
genus Betacoronavirus. 

Nucleotide sequences of partial RdRp-helicase genes 
derived from fecal samples IFB2012-17F and IFB2012-8F 
were deposited under accession numbers AB918718 and 
AB918719, respectively, in the DNA Data Bank of Japan 
(DDBJ, Tokyo) nucleotide database. 

Various tissues and tracheal swab samples collected 
from the bats where the feces yielded a positive result were 
also subjected to nested RT-PCR for the partial RdRp 
region. Spleen tissue of IFB2012-13F, brain tissue of 
IFB2012-17F, and lung tissue of IFB2012-17F showed 
positive signals for RdRp gene amplification. Sequence 
analysis of tissue samples showed sequences identical to 
the partial RdRp sequence obtained from fecal samples 
(IFB2012-13F and IFB2012-17F). The respiratory epithe- 
lium is known to be a target for CoVs [27], and therefore, 
detection of bat CoV genome fragments in lung samples 
might be expected. However, the detection of CoV genome 
fragments in spleen and brain samples suggested that bat 
CoV infection in these species may not be restricted to the 
respiratory or alimentary tract. 

Attempts at virus isolation using Vero E6, FBKT, and 
DemKT1 cells were unsuccessful. No CPE was observed in 
cell cultures, and viral RNA was not detected by RT-PCR 
from the culture supernatants after the fifth passage. 

Nucleotide sequencing and phylogenetic analysis of 
partial RdRp-helicase regions of the samples in the present 
study and these of known CoVs suggest that the CoV 
genome fragments obtained in this study are related to 


Location No. of samples No. of positive 
tested samples 
Paguyaman 17 3 
18 0 
2 0 
Surabaya 3 0 
Yogyakarta 19 0 
Paguyaman 7 0 
8 0 
74 3 
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Fig. 2 Phylogenetic tree of partial RdRp-helicase using a 2,116-bp 
nucleotide sequence following removal of primer sequences at the 5’ 
and 3’ end of the partial RdRp-helicase gene obtained from fruit bats 
samples (IFB2012-8F and IFB2012-17F) and various coronaviruses 
derived from GenBank. The tree was constructed by the maximum- 


those of members of the genus Betacoronavirus. More than 
90 % amino acid sequence identity in the conserved rep- 
licase domain, ADRP, nsp5 (3CLpro), nsp12 (RdRp), 
nsp13 (helicase), nspl4 (ExoN), nsp15 (NendoU), and 
nsp16 (O-MT) should be observed in order to include CoVs 
in the same species [28]. Therefore, the CoVs detected in 
this study could not be defined at the species level, since 
only certain regions of conserved replicase domains were 
determined. However, the number of nucleotide substitu- 
tions over the partial RdRp-helicase sequence between 
IFB2012-8F and IFB2012-17F and phylogenetic analysis 
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likelihood method using MEGA 6 software. Bootstrap values were 
calculated for 1000 replicates and are shown next to the branches. Bat 
coronaviruses detected in this study are indicated by black squares. 
SARS-CoV, severe acute respiratory syndrome coronavirus; avian 
IBV, avian infectious bronchitis virus; CoV, coronavirus 


of partial RdRp-helicase suggested that the fruit bats harbor 
variants of a single coronavirus strain. 

In this study, CoV genomes containing regions resem- 
bling those of members of the genus Betacoronavirus were 
detected in Dobsonia moluccensis, known as Moluccan 
naked-backed fruit bat, in Indonesia. Dobsonia moluccen- 
sis belongs to family Pteropodidae, which is known to 
harbor a number of viruses, including CoVs [25, 29], pa- 
ramyxoviruses [18, 30] and rhabdoviruses [29, 31]. Nipah 
and Hendra viruses have also been detected in Dobsonia 
moluccensis [32, 33]. 


Coronavirus genomes in fruit bats in Indonesia 


This report highlights the fact that Indonesian fruit bats 


are a natural reservoir for various viruses, which have the 
potential to be zoonotically transmitted, since fruit bats are 
also consumed in certain parts of Indonesia. 
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